Carrier Density Dependence of the Magnetic Properties in Iron-doped Bi2Se3 
Topological Insulator 

H. Li,-^ Y. R. Song^ Meng-Yu Yao,^ Fengfeng Zhu/ Canhua Liu,^ C. L. Gao,^ Jin-Feng 
Jia,^ Dong Qian^y X. Yao.^ Y. J. Shi,^ and D.Wu^ 

^^Key Laboratory of Artificial Structures and Quantum Control (Ministry 
of Education), Department of Physics, Shanghai Jiao Tong University, 
Shanghai 200240, China 

"^^ The National Laboratory of Solid State Micro structures (NLSSM), 
Nanjing University, Nanjing, 210093, China 

The electronic and magnetic properties of iron-doped topological insulator 
Bii.84-xFeo.i6Ca2:Se3 single crystals were studied. By co-doping Fe and Ca atoms, 
ferromagnetic bulk states with different carrier density (from n-type to p-type) were 
obtained. Effective magnetic moments for each Fe atom was estimated as small as 
about 0.07//^. Magnetic and non-magnetic phase separation was observed in all 
samples. Our results suggest that the bulk ferromagnetism in Fe-doped Bi2Se3 is not 
intrinsic and regardless of carrier density. 
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I. INTRODUCTION 



Distinguished from trivial insulators, topological insulators (TIs) - a new state of quan- 
tum matter - possessing gapless topological edge states protected by time reversal symmetry 
(TRS) have been attracting extensive studies recently^"— . A number of materials have been 
proved or proposed to be TIs including Hg(Cd)Te quantum well^, Bii.a-Sbj; alloys^, binary 
compounds (Bi2Se3, Bi2Te3, Sb2Te3}^""— i^^, Half-Heusler compounds^^i^, ultra-thin Bi(lll) 
bilayersi^ii^ and so on. Among these materials, Bi2Se3 is thought to be the most promis- 
ing three dimensional TI for potential applications because of its simplest single helical 
Dirac cone surface states and the largest bulk gap (300meV)^i^. An energy gap can be 
opened at Dirac point to form massive Dirac Fermion by introducing TRS breaking field, 
such as magnetic field. Similar to graphene, opening an energy gap in the Dirac cone is 
very important for potential applications. In addition, the interplay between ferromag- 
netic states and topological states can cause many exotic phenomenon such as magnetic 
monopole, anomalous quantum Hall effect and the topological properties related Faraday 
and Kerr magneto-optical effects^""—. Ferromagnetic states have been observed in 3d tran- 
sitional magnetic elements (V, Cr, Mn, Fe) doped Bi2Te3, Bi2Se3, Sb2Te3 single crystals or 
filma^"— . Massive Dirac Fermion will form when the magnetic moments have component 
perpendicular to the surface. In those magnetic TIs, the expected massive Dirac Fermion 
was only observed in Fe-doped and Mn-doped Bi2Se3 using angle-resolved photoemission 
spectroscopy (ARPES)^i^. Though, the observed energy gap at Dirac point is very large 
(~60meV), the bulk magnetic properties of Fe^Bi2_xSe3 are still controversial2^i22i2^. At the 
same time, the Fermi level is not within the gap^ either. In the first report which claimed 
the observation of the massive Dirac Fermion^, high Curie temperature (Tc) and small 
total magnetic moments were observed. Assuming doped Fe atoms have a valence of +3, 
only few Fe (x=0.002 - 0.005 ) atoms are estimated to enter the material matrix with a 
nominal concentration of x=0.1 to x=0.3. Intergrowth of Bi2Se3 and magnetic Fe^Se^^ phase 
was also observed in some samples^. In another experiment^, Fe atoms were thought to 
completely enter the material matrix and make Fe-doped Bi2Se3 be intrinsic ferromagnetic 
semiconductor. Since Bi2Se3 is a semiconductor, in principle, intrinsic ferromagnetic TI 
would have similar magnetic behavior as the well-defined diluted magnetic semiconductor 
such as (Mn,Ga)As2ii2^. If the ferromagnetic states in Fe-doped Bi2Se3 is intrinsic, the over- 
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all magnetic properties should strongly depend on the carrier density. On the other hand, 
in some systems, dopants with different valance can destroy the formation of ferromagnetic 
clusters^^. It is very important to carry out the carrier dependent experiments on Fe-doped 
Bi2Se3 systems for several reasons: i)To understand the origin of the bulk ferromagnetic 
states. ii)To control its magnetic properties. iii)To move the Fermi level into the gap near 
the Dirac point. 

In this work, we carefully studied the magnetic properties of Fe-doped Bi2Se3 with Fe 
concentration of x=0.16 where the energy gap at Dirac point was nearly saturated^- as a 
function of carrier density. We succeeded in tuning the Fermi level to Dirac point with 
proper mount of Ca doping. For all Ca doped samples, ferromagnetic states were obtained. 
However, magnetic and non-magnetic phases separation existed in all samples. The magnetic 
properties could not be tuned systematically by carrier density. Regardless the carrier 
density, bulk ferromagnetism mainly came from extrinsic magnetic clusters. 

II. EXPERIMENT 

High quality single crystals of Bii.84-zFeo.i6Ca2,.Se3 (x=0, 0.02, 0.04, 0.06) were grown 
by modified Bridgman method^. High purity of Bi (99.999%), Fe (99.999%), Ca (99.5%) 
and Se (99.999%) powders were carefully mixed and sealed in evacuated quartz tubes. Ex- 
cess amount of selenium was used to reduce the selenium vacancy. The tubes were firstly 
heated to 850°C and kept for 24 hours, then slowly cooled to 550 °C within 2 days, fol- 
lowed by 3 days annealing. The obtained crystals were well crystallized and could be easily 
cleaved perpendicular to the c-axis. The concentration of Fe and Ca dopants in the crys- 
tals were carefully determined by inductively coupled plasma (ICP) (iCAP THERMO 6000 
Radial). The band structures were measured by ARPES with 10 eV photons at 77K in 
National Synchrotron Radiation Laboratory (NSRL, Hefei) using Scienta R4000 analyzer 
with base pressures better than 5x10^^^ Torr. Hall measurements were carried out in the 
Magnetocryostat (Oxford). Magnetic measurements were performed by Physical Property 
Measurement System (PPMS) (Quantum Design) with magnetic field paralleled or perpen- 
dicular to the c-axis of the samples. ICP analysis shows that almost all of the nominal 
Fe atoms incorporated into the crystals. Though not all the Ca can enter the crystals, Ca 
concentration increases with x, approximately (Table-1). The maximum Ca entered into the 
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crystals we obtained is about Xe//=0.018. 



III. RESULTS AND DISCUSSION 

As-grown Bi2Se3 crystals show electron doping (n-type) due to the intrinsic Se vacancies^i^. 
Altering the Bi:Se ratio in the nominal composition can reduce the defects to a certain ex- 
tent, however, can not change the carrier type. If Fe atoms enter the Bi-Se matrix, they can 
occupy two possible positions in the crystal. The first possibility is that Fe atoms substitute 
Bi atoms with the same valence of +3, thus it will hardly change the carrier density^i^i^. 
This kind of substitution happens in other doped topological insulators2^i2£i2£. Another 
possibility is Fe intercalation in the Van der Waals space between two quintuple layers. In 
this case, Fe atoms will supply additional electrons into Fermi seai^i^. Figure 1 (a)-(d) 
present the low energy electronic structures of Bii.84-xFeo.i6Ca2:Se3 measured by ARPES at 
77K. Linearly dispersive surface states are clearly observed in all samples. Similar to the 
reported Fe-doped samplesSS, the Dirac point of surface states locates at ~ 200 meV below 
the Fermi level without Ca doping. Fermi level touches the bottom of conduction bands. In 
Bi2Se3 single crystals, the Dirac point locates at about 200-350 meV below the Fermi level 
depending on the initial ratio of Bi and Se. Known from the ARPES results, Fe does not act 
as an electron donor in our samples, which indicates that there is very little Fe intercalation. 
Each Ca atom can substitute Bi atom in the form of Ca^+ and supplies a hole into the Fermi 
seaF^, so with the increase of Ca doping levels the Fermi level moves downward. Fermi level 
is tuned to the position that is very close the Dirac point (Fig. 1(d)) when the Ca nominal 
doping is x=0.06. Figure 1 (e) and (f) show the energy distribution curves (EDCs) near the 
Dirac points for x=0 and x=0.02, respectively. Blue curves are the EDCs at k=0. Spectral 
weight suppression ("gap"-like) at Dirac point were clearly observed on x=0 and x=0.02 
samples. Seen from Fig. 1(e) and (f), though the energy resolution is not as good as the 
reference [30], two peaks in EDC (blue curve) were observed which is very similar to the 
previous report^. For x=0.04 and 0.06 samples, because of the doping induced disorder, the 
spectra is not sharp enough to observe spectral weight suppression. The gap we observed in 
x=0 and x=0.02 samples is about 50 to 90 meV depending on samples that is close to the 
gap value in reference [30]. However the "gap"-like feature is not the sole evidence of the 
formation of long-range intrinsic ferromagnetism. In fact, "gap" -like feature was observed 
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FIG. 1. Low energy ARPES spectra of Bii.84--xFeo.i6Ca2,Se3 with (a)x=0, (b)x=0.02, (c)x=0.04, 

(d) x=0.06. Fermi level moves downwards to the Dirac point with the increase of Ca dopants. 

(e) EDCs near the Dirac point for x=0 and (f)x=0.02. Blue curves are EDCs at k=0. (g) Hall 
resistivity as a function of applied magnetic filed, (h) Carrier density as a function of Ca concen- 
tration. 

in many topological insulators with magnetic or non-magnetic doping recently—. We can 
not make conclusion on the magnetic states based on the ARPES results. ARPES mea- 
surement is surface sensitive, which can give the information of the carrier density near the 
surface region. The bulk carrier density as a function of Ca doping was further measured 
by Hall effects. Figure 1(g) presents the Hall resistivity as a function of applied magnetic 
field. Figure 1(h) shows the carrier density extracted from Fig. 1(g) as a function of Ca 
concentration. The reported carrier density of Bi2Se3 varies from ~10^^ cm^^ to ~10^'' cm~^ 
obtained by different groups^"— depending on the initial stoichiometry of Bi and Se and the 
preparing methods. In our parent compound Bii 84Feo.i6Se3, the measured carrier density 
is about 2.5xl0^^cm~^ which is consistent with previous reports. Consistent with ARPES 
measurements, the incorporation of Ca into the crystals reduces the electron density and 
drives the carrier change from n-type to p-type. By co-doping Fe and Ca, we successfully 
tuned the Fermi level to the Dirac point. The successful tuning of the carrier concentration 
also allows us to explore the relation between magnetic properties and carrier density. 

For all the samples, bulk ferromagnetism is observed at 2. IK. Figure 2 (a) and (b) show 
the magnetic moments (M) as a function of temperature (T) for different Ca concentrations 
with magnetic field perpendicular or parallel to the c-axis, respectively. Under the applied 
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FIG. 2. Field cooled temperature dependent dc magnetism with an applied filed of IkOe for the 
Bii.84-a;Feo.i6Caa;Se3 crystals along (a) in-plane, (b)out-of-plane, (c)in-plane after normalization, 
(d) out-of-plane after normalization. 



magnetic filed (IKOe), the amplitude of the magnetic signal varies for different Ca concen- 
tration, however, they show similar T-dependent behavior. Figure 2 (c) and (d) present 
the M-T curves after being normalized by the signals at the lowest temperature (2.1 K). 
In other magnetic topological insulators, like Cr-doped Sb2Te3^'2E and Mn-doped Bi2Te^, 
only single magnetic transition point was observed. However, in Fe-doped Bi2Se3, there are 
three possible turn points in the M-T curves - at about lOK, lOOK to 200K and 250K to 
300K, similar to previous results^'^. This implies the magnetic properties in this system 
are much more complicated. 

In order to find the origin of the magnetic signals in Fe-doped Bi2Se3, we further sys- 
tematically studied the magnetic hysteresis loops (M-H curves) along the in-plane direction 
as a function of carrier density. Figure 3 present M-H curves of samples with different Ca 
concentrations. The M-H curve of Bi2Se3 sample shows pure diamagnetic behavior, which 
does not show any noticeable change by small amount of Ca-doping. The blue curves present 
diamagnetic signals from the parent compounds without Fe. With Fe doping, there are clear 
hysteresis loops at 2. IK which indicates the existence of ferromagnetic state. Seen from Fig. 
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FIG. 3. The Magnetic hysteresis loops along the in-plane direction at different temperature for 
(a)x=0, (b)x=0.02 (c)x=0.04, (d)x=0.06 
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FIG. 4. The Magnetic hysteresis loops along the in-plane direction at different temperatures after 
subtracting the diamagnetic signals for (a)x=0, (b)x=0.02 (c)x=0.04, (d)x=0.06 
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FIG. 5. Comparison of the magnetic hysteresis loops along the out-of-plane and in-plan directions 
after subtracting the diamagnetic signals for (a)x=0, (b)x=0.02 (c)x=0.04, (d)x=0.06 

3(a), the magnetic signal is nearly saturated at about 40 KOe and a diamagnetic signal 
appears again at high magnetic field. Insert in Fig. 3(a) presents the magnetic signal up 
to 90 KOe. Diamagnetic signals are much clearer at high field. The estimated diamagnetic 
susceptibility is about 7.1xlO~®emu/(g-T) which is close to that of pure Bi2Se3. Similar 
diamagnetic signals are also obvious in Ca-doped samples (Fig. 3(b)-(d)). The existence of 
the diamagnetic signal implies that there is phase separation of non-magnetic Bi2Se3 or Ca- 
doped Bi2Se3 phase and magnetic phase. In Figure 4, We subtract the diamagnetic signals 
from the total signals to get the pure magnetic signals that related to the magnetic phase 
with Fe dopants. The overall M-H curves become simple. See from Fig. 5, at 2. IK, at zero 
field, the remanent magnetization is larger along the in-plane direction, which means the 
magnetic moments intend to stay in plane. 

For all the samples, the pure magnetic signals are nearly saturated at 60 KOe. Using the 
signal measured at 60 KOe and ICP results, we calculated the average magnetic moments of 
per Fe atom as shown in Table 1. For all samples, the calculated magnetic moments per Fe 
atom are very close. It is about 0.07±0.02y[XB which is much smaller than possible magnetic 
moments of Fe^"^ ion. Such small magnetic moment can not been explained by the scenario 
of Fe substitution of Bi. Actually, first principle calculation^^ assuming Fe substitution of 
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TABLE I. Concentrations of doped Ca and Fe and the magnetic moments for different magnetic 
dopants 



X 


Fe 


Ca 


n(cm ^) 


M-QTifJ-B/Featom) 





0.16±0.02 





-2.47x1018 


O.OSibO.Ol 


0.02 


0.19±0.02 


0.009±0.005 


-7.04x10^^ 


0.07±0.01 


0.04 


0.16±0.02 


0.012±0.005 


-4.46 xlO^^ 


0.06±0.01 


0.06 


0.14±0.02 


0.018±0.005 


2.14xl0i'^ 


0.06±0.01 



Bi gave high spin states of Fe atoms. Since the magnetic properties of Fe-doped Bi2Se3 do 
not obviously depend on the carrier density, we conclude that the most part of the magnetic 
signals is extrinsic. They are not diluted ferromagnetic semiconductors. We can exclude the 
formation of Fe inclusions. If the ferromagnetism comes from Fe inclusions, we should get 
large average magnetic moments per Fe atoms(2.2/xs/atom for bulk Fe). On the other hand, 
Fe and Se can form complicated compounds Fe^Sey—^ which present ferrimagnetic and/or 
antiferromagnetic behaviors. According to previous report on bulk Fe^^Scy compounds^, 
the value of magnetic moment of Fe atom in Fe^^Scy varies depending on the exactly ration 
of Fe and Se. The reported minimum value is from 0.2 to 0.4 fis and has the same order 
of magnitude as our results. So we think that the magnetic phase in our Fe-doped Bi2Se3 
samples could be some forms of FejjSej^ nanocrystals. Two possible explanation on the very 
small magnetic momentum. First, those nanocrystals are ferrimagnetic with low spin states. 
Second, they can be in antiferromagnetic states with uncompensated surface spins that gives 
the small ferromagnetic signals. Since our finding is based on bulk measurement, we can not 
exclude some intrinsic ferromagnetism on the surface which may open a real gap in ARPES 
spectra. 

IV. CONCLUSIONS 

By co-doping Ca and Fe, we succeeded in tuning the carrier density from n-type to p-type 
in ferromagnetic Bii,84-xFeo.i6Caa;Se3 samples. We studied the evolution of the magnetic 
properties as a function of carrier density. Contrast to well-defined ferromagnetic semicon- 
ductor, the magnetic properties in Fe-doped Bi2Se3 do not have obvious dependence on the 
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carrier density. Our results suggest that the bulk ferromagnetic behaviors of Fe-doped Bi2Se3 
are due to FeajSe^ nanoclusters in the crystals. In principle, the nature of the gap observed 
near Dirac point could be solved by determining the spin polarization using spin-resolved 
ARPES or spin-polarized scanning tunnelling spectroscopy in the future. 
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